CMB Bispectrum from Active Models of Structure Formation 
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We propose a new method for a numerical computation of the angular bispectrum of the CMB 
anisotropies arising from active models such as cosmic topological defects, using a modified Boltz- 
mann code based on CMBFAST. The method does not use CMB sky maps and requires moderate 
computational power. As a first implementation, we apply our method to a recently proposed model 
of simulated cosmic strings and estimate the observability of the non-Gaussian bispectrum signal. 
A comparison with the cosmic variance of the bispectrum estimator shows that the bispectrum for 
the simulated string model we used is not observable. 



Anisotropies of the Cosmic Microwave Background ra- 
diation (CMB) are directly related to the origin of struc- 
ture in the universe. Galaxies and clusters of galaxies 
eventually formed by gravitational instability from pri- 
mordial density fluctuations, and these same fluctuations 
left their imprint on the CMB. Recent balloon |j]f| and 
ground-based interferometer experiments have pro- 
duced reliable estimates of the power spectrum of the 
CMB temperature anisotropies. While they helped elim- 
inate certain candidate theories for the primary source 
of cosmic perturbations, the power spectrum data is still 
compatible with the theoretical estimates of a relatively 
large variety of models, such as ACDM, quintessence 
models or some hybrid models including cosmic defects. 
These models, however, differ in their predictions for 
the statistical distribution of the anisotropies beyond the 
power spectrum. Future MAP and Planck satellite mis- 
sions (scheduled for launch in 2001 and 2007, respec- 
tively) will provide high-precision data allowing definite 
estimates of non- Gaussian signals in the CMB. It is there- 
fore important to know precisely which are the predic- 
tions of all candidate models for the statistical quantities 
that will be extracted from the new data and identify 
their specific signatures. 

There are two main classes of models of structure 
formation — passive and active models. In passive mod- 
els, density inhomogencitics are set as initial conditions 
at some early time, and while they subsequently evolve 
as described by Einstein-Boltzmann equations, no addi- 
tional perturbations are seeded. On the other hand, in 
active models the sources of density perturbations are 
time-dependent . 

All specific realizations of passive models are based 
on the idea of inflation. In most inflationary models, 
density fluctuations arise from quantum fluctuations of 
a scalar field placed in the vacuum and hence are well 
described by a Gaussian distribution, while second-order 
effects may add a weak non-Gaussian signal [||. 
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On the other hand, active models of structure forma- 
tion are motivated by cosmic topological defects. If our 
ideas about grand unification are correct, then some cos- 
mic defects, such as domain walls, strings, monopoles or 
textures, should have formed during phase transitions in 
the early universe 0. Once formed, cosmic strings could 
survive long enough to seed density perturbations. De- 
fect models possess the attractive feature that they have 
no parameter freedom, as all the necessary information is 
in principle contained in the underlying particle physics 
model. Generically, perturbations produced by active 
models are not expected to be Gaussian distributed. 

The narrow main peak and the presence of the second 
and the third peaks in the CMB angular power spec- 
trum, as measured by BOOMERANG, MAXIMA and 
DASI , is an evidence for coherent oscillations of the 
photon-baryon fluid at the beginning of the decoupling 
epoch JtJ . While such coherence is a property of all pas- 
sive model, realistic cosmic string models produce highly 
incoherent perturbations that result in a much broader 
main peak. This excludes cosmic strings as the primary 
source of density fluctuations unless new physics is pos- 
tulated, e.g. varying speed of light B. In addition to 
purely active or passive models, it has been recently sug- 
gested that perturbations could be seeded by some com- 
bination of the two mechanisms. For example, cosmic 
strings could have formed just before the end of infla- 
tion and partially contributed to seeding density fluctua- 
tions. It has been shown || that such hybrid models can 
be rather successful in fitting the CMB power spectrum 
data. Therefore, statistics beyond the power spectrum 
is required to discriminate between active and passive 
models. 

Of the available non-Gaussian statistics, the CMB bis- 
pectrum, or the three-point function of Fourier compo- 
nents of the temperature anisotropy, has been perhaps 
the one best studied in the literature pO|Jll]] . Although 
there are a few cases where the bispectrum may be esti- 
mated analytically from the underlying model, a precise 
numerical code to compute it, similar to the CMBFAST 
code [Q for the power spectrum, is presently lacking. 
The bispectrum can be estimated from simulated CMB 
sky maps; however, computing a large number of full-sky 
maps resulting from defects is a much more demanding 
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task. 

In this article we introduce a new method for obtaining 
the CMB bispectrum directly from numerically simulated 
defect models, without building CMB sky maps. Given 
a suitable model, one can generate a statistical ensemble 
of realizations of defect matter perturbations. We use a 
modified Boltzmann code based on CMBFAST to com- 
pute the effect of these perturbations on the CMB and 
find the bispectrum estimator for a given realization of 
sources. We then perform statistical averaging over the 
ensemble of realizations to compute the expected CMB 
bispectrum. (The CMB power spectrum is also obtained 
as a byproduct.) Our method is specifically tailored for 
computations of the bispectrum; extending it to higher- 
order correlation functions would require prohibitively 
longer calculations. As a first application, we computed 
the expected CMB bispectrum from a model of simulated 
string networks first introduced by Albrecht et al. C3] 
and further developed in Ref. jl4| and in this work. Our 
calculations indicate that the bispectrum resulting from 
this model is negligible when compared with the cosmic 
variance. We discuss the implications of this result for 
detectability of cosmic strings through the bispectrum 
statistic. 



I. CMB BISPECTRUM FROM ACTIVE MODELS 

We assume that, given a model of active perturbations, 
such as a string simulation, we can calculate the energy- 
momentum tensor T^^x, r) for a particular realization 
of the sources in a finite spatial volume Vq. Here, x is a 
3-dimensional coordinate and r is the cosmic time. Many 
simulations are run to obtain an ensemble of random real- 
izations of sources with statistical properties appropriate 
for the given model. The spatial Fourier decomposition 
of can be written as 

T^(x,r)=^e Al ,(k,r)e jkx , (1) 

k 

where k are discrete. If Vq is sufficiently large we can 
approximate the summation by the integral 

£ e^(k, r)e* kx « ^ J d 3 kO^(k, r)e 4kx , (2) 

and the corresponding inverse Fourier transform will be 
e^(k,r) = i/ d 3 xr F (x,r) e - kx . (3) 

V JVo 

Of course, the final results, such as the CMB power spec- 
trum or bispectrum, do not depend on the choice of Vq. 
To ensure this independence, we shall keep Vq in all ex- 
pressions where it appears throughout the following sec- 
tions. 

It is conventional to expand the temperature fluctua- 
tions over the basis of spherical harmonics, 



AT/T(n) = ^a im F /m (n), (4) 

Irn 

where n is a unit vector. The coefficients ai m can be 
decomposed into Fourier modes, 

aim = jBp 477 J d " kAi (k) YL[k) - (5) 

Given the sources jUI/ (k, t), the quantities A;(k) are 
found by solving linearized Einstein-Boltzmann equa- 
tions and integrating along the line of sight, using a 
code similar to CMBFAST @. This standard proce- 
dure can be written symbolically as the action of a linear 
operator B^ u (k) on the source energy-momentum ten- 
sor, A;(k) = B^ u {k)Q ^{k, r), so the third moment of 
A;(k) is linearly related to the three-point correlator of 
(k, r) . Below we consider the quantities A; (k) , cor- 
responding to a set of realizations of active sources, as 
given. The numerical procedure for computing A;(k) 
was developed in Refs. [jl3| and Q. 

The third moment of a im , namely (ai irni ai 2m2 ai 3m3 ) , 
can be expressed as 

(-i) 1 ^ 1 " (4 7 r) 3 ^| ( i 3 k 1 d 3 k 2( i 3 k3Y^ mi (k 1 ) 

x Y*: m2 (k 2 )Y£ m3 (k 3 ) (A^kx) A i2 (k 2 ) A i3 (k 3 )} . (6) 

A straightforward numerical evaluation of Eq. (||) from 
given sources A; (k) is prohibitively difficult, because it 
involves too many integrations of oscillating functions. 
However, we shall be able to reduce the computation to 
integrations over scalars (a similar method was employed 
in p"5[). Due to homogeneity, the 3-point function van- 
ishes unless the triangle constraint is satisfied, 

ki + k 2 + k 3 = 0. (7) 

We may write 

(A h (k x )A i2 (k 2 )A i3 (ks)> 
= <5 (3) (k x + k 2 + ka) P hhh (ki,k 2 , k 3 ) , (8) 

where the three-point function Pi 1 i 2 i 3 (ki,k 2 ,k3) is de- 
fined only for values of k^ that satisfy Eq. (^) . Given the 
scalar values k\ , fc 2 , fc 3 , there is a unique (up to an over- 
all rotation) triplet of directions k^ for which the RHS of 
Eq. (||) does not vanish. The quantity P/ 1 ; 2 / 3 (ki,k 2 ,k 3 ) 
is invariant under an overall rotation of all three vectors 
ki and therefore may be equivalcntly represented by a 
function of scalar values k\, fc 2 , & 3 , while preserving all 
angular information. Hence, we can rewrite Eq. (||) as 

(A ;i (k 1 )A i2 (k 2 )A ;3 (k 3 )) 
= <S (3) (ki + k 2 + k 3 ) P h i 2h {k t , k 2 , ki). (9) 

Then, using the simulation volume Vq explicitly, we have 
Pi Ms (*i,fca, fcs) = (A^kx) A i2 (k 2 ) A i3 (k 3 )) . (10) 
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Given an arbitrary direction ki and the magnitudes fci, 
k 2 and k 3 , the directions k 2 and k3 are specified up to 
overall rotations by the triangle constraint. Therefore, 
both sides of Eq. (|Io| ) are functions of scalar fcj only. 
The expression on the RHS of Eq. ( [To| ) is evaluated nu- 
merically by averaging over different realizations of the 
sources and over permissible directions kj ; below we shall 
give more details of the procedure. 

Substituting Eqs. (^) and ([lO]) into (^), Fourier trans- 
forming the Dirac delta and using the Rayleigh identity, 
we can perform all angular integrations analytically and 
obtain a compact form for the third moment, 



Kffi"* I r'drb hl2l3 (r), (11) 



h h h 



where, denoting the Wigner 3j-symbol by V)l) in i(iJ . 
have 
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and where we have defined the auxiliary quantities bi t i 2 i 3 
using spherical Bessel functions ji, 



hil 2 l 3 (r) 
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7T 3 (27r) 3 



k\dk\ k\dk 2 k\dk 3 



x jh{hr)ji 2 (k2r)ji 3 {k 3 r)Pi 1 i 2 i 3 {ki,k 2 ,k 3 ). (13) 



The volume factor V contained in this expression is 
correct: as shown in the next section, each term A/ 



includes a factor V 2 ^ 3 , while the average quantity 
Pi 1 i 2 i 3 (ki,k 2 ,k 3 ) oc V ~ 3 [cf. Eq. @], so that the ar- 
bitrary volume Vb of the simulation cancels. 

Our proposed numerical procedure therefore consists 
of computing the RHS of Eq. (|lTJ) by evaluating the nec- 
essary integrals. For fixed {hl 2 l 3 }, computation of the 
quantities buuu(r) is a triple integral over scalar fcj de- 
fined by Eq . (|13| ) ; it is followed by a fourth scalar integral 
over r [Eq. flllD]. We also need to average over many real- 
izations of sources to obtain Pi x i 2 i 3 (fci, k 2 , k 3 ). It was not 
feasible for us to precompute the values P^i^ (fci, k 2 , k 3 ) 
on a grid before integration because of the large vol- 
ume of data: for each set {^1^2/3} the grid must contain 
~ 10 3 points for each fcj. Instead, we precompute A;(k) 
from one realization of sources and evaluate the RHS of 
Eq. @ on that data as an estimator of Pi t i 2 i 3 (fci, k 2 , k 3 ), 
averaging over allowed directions of kj . The result is used 
for integration in Eq. (13). 



Because of isotropy and since the allowed sets of di- 
rections ki are planar, it is enough to restrict the nu- 
merical calculation to directions k; within a fixed two- 
dimensional plane. This significantly reduces the amount 
of computations and data storage, since A;(k) only needs 
to be stored on a two-dimensional grid of k. 



In estimating Pi 1 i 2 i 3 (kx, A^, fca) from Eq. fllO|), averag- 
ing over directions of k, plays a similar role to ensemble 
averaging over source realizations. Therefore if the num- 
ber of directions is large enough (we used 720 for cosmic 
strings), only a moderate number of different source re- 
alizations is needed. The main numerical difficulty is the 
highly oscillating nature of the function bi x i 2 i 3 (r). The 
calculation of the bispectrum for cosmic strings presented 
in Section [H] requires about 20 days of a single-CPU work- 
station time per realization. 

We note that this method is specific for the bispectrum 
and cannot be applied to compute higher-order correla- 
tions. The reason is that higher-order correlations in- 
volve configurations of vectors k^ that are not described 
by scalar values ki and not restricted to a plane. For 
instance, a computation of a 4-point function would in- 
volve integration of highly oscillating functions over four 
vectors k^ which is computationally infeasible. 

From Eq. ( [ji"] ) we derive the CMB angular bispectrum 
C h i 2 i 3 , defined as 0] 



(ai imi ai 



2fll2 <^3m3 



h h h 
mi m 2 m 3 



hhh 



(14) 



The presence of the 3j-symbol guarantees that the third 
moment vanishes unless mi + m 2 + m 3 = and the U 
indices satisfy the triangle rule \h ~ lj \ < h < h + lj- In- 
variance under spatial inversions of the three-point cor- 
relation function implies the additional 'selection rule' 
li + l 2 + l 3 = even, in order for the third moment not to 
vanish. Finally, from this last relation and using standard 
properties of the 3j-symbols, it follows that the angular 
bispectrum Ci t i 2 i 3 is left unchanged under any arbitrary 
permutation of the indices Zj. 

In this paper we restrict our calculations to the angular 
bispectrum Ci^i 3 in the 'diagonal' case, i.e. l\ = l 2 = 
l 3 = I. This is a representative case and, in fact, the one 
most frequently considered in the literature. Plots of the 
power spectrum are usually done in terms of 1(1 + 1)C; 
which, apart from constant factors, is the contribution to 
the mean squared anisotropy of temperature fluctuations 
per unit logarithmic interval of I. In full analogy with 
this, the relevant quantity to work with in the case of 
the bispectrum is 



G 



HI 



l(2l + l) 3/2 



C, 



in 



(15) 



Details of this derivation are presented in the Appendix. 
For large values of the multipole index I, Gm oc l 3 ^ 2 Cm. 
Note also what happens with the 3j-symbols appearing 



in the definition of the coefficients TC 



hhh 



: the symbol 



( h h h \ ; 

\m-\ 7772771 3 / 

in Eq. © 



is absent from the definition of C, 



while 



the symbol (000) is squared. Hence, there 

are no remnant oscillations due to the alternating sign of 

(i 11 \ 

\oooj- 

However, even more important than the value of Cm it- 
self is the relation between the bispectrum and the cosmic 
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variance associated with it. In fact, it is their comparison 
that tells us about the observability 'in principle' of the 
non-Gaussian signal. The cosmic variance constitutes a 
theoretical uncertainty for all observable quantities and 
comes about due to the fact of having just one realization 
of the stochastic process, in our case, the CMB sky \ ^T\ . 

The way to proceed is to employ an estimator Ci ± i 2 i a 
for the bispectrum and compute the variance from it. 
By choosing an unbiased estimator we ensure it satisfies 
QiWa = (QiWs)- However, this condition does not iso- 
late a unique estimator. The proper way to select the best 
unbiased estimator is to compute the variances of all can- 
didates and choose the one with the smallest value. The 
estimator with this property was computed in Ref. |lq] 
and is 

Chhh = E ( m x m 2 m 3 ) a hm 1 ai 2 m 2 ai ama . (16) 

mi ,7712 ,7713 ^ ' 

The variance of this estimator, assuming a mildly non- 
Gaussian distribution, can be expressed in terms of the 
angular power spectrum Ci as follows jll| 

0> = 0^012012(1 + 8^1.,+ Si 2 i 3 + Si 3 i 1 +25i 1 i 2 Si 2 i 3 ) . (17) 

The theoretical signal-to-noise ratio for the bispectrum 
is then given by 

(s/N) hl2l3 = \c hhl3 /a dHi2i j, (is) 

In turn, for the diagonal case l± — I2 — h — I we have 

(S/N)i = \C m /a dm \. (19) 

Incorporating all the specifics of the particular exper- 
iment, such as sky coverage, angular resolution, etc., 
will allow us to give an estimate of the particular non- 
Gaussian signature associated with a given active source 
and, if observable, indicate the appropriate range of mul- 
tipole Us where it is best to look for it. 

II. BISPECTRUM FROM STRINGS 
A. The string model 

To calculate the sources of perturbations we use an up- 
dated version of the cosmic string model first introduced 
by Albrecht et al. JjlJ and further developed in Ref. Q , 
where the wiggly nature of strings was taken into ac- 
count. In these previous works the model was tailored to 
the computation of the two-point statistics (matter and 
CMB power spectra). When dealing with higher-order 
statistics, such as the bispectrum, a different strategy 
needs to be employed. 

In the model, the string network is represented by a 
collection of uncorrelated straight string segments pro- 
duced at some early epoch and moving with random un- 
correlated velocities. At every subsequent epoch, a cer- 
tain fraction of the number of segments decays in a way 



that maintains network scaling. The length of each seg- 
ment at any time is taken to be equal to the correlation 
length of the network. This and the root mean square 
velocity of segments are computed from the velocity- 
dependent one-scale model of Martins and Shellard [jl8| . 
The positions of segments are drawn from a uniform dis- 
tribution in space, and their orientations are chosen from 
a uniform distribution on a two-sphere. 

The total energy of the string network in a volume V 
at any time is E = Np,L, where N is the total num- 
ber of string segments at that time, [i is the mass per 
unit length, and L is the length of one segment. If 
L is the correlation length of the string network then, 
according to the one-scale model, the energy density is 
p = E/V = f-i/L 2 , where V = Vba 3 , the expansion factor 
a is normalized so that a = 1 today, and Vq is a constant 
simulation volume. It follows that N = V/L 3 = Vo/£ 3 , 
where £ = L/a is the comoving correlation length. In the 
scaling regime £ is approximately proportional to the con- 
formal time r and so the number of strings N(t) within 
the simulation volume Vo falls as r~ 3 . 

To calculate the CMB anisotropy one needs to evolve 
the string network over at least four orders of magnitude 
in cosmic expansion. Hence, one would have to start 
with N > 10 12 string segments in order to have one seg- 
ment left at the present time. Keeping track of such a 
huge number of segments is numerically infeasible. A way 
around this difficulty was suggested in Ref. |^3| , where 
the idea was to consolidate all string segments that decay 
at the same epoch. The number of segments that decay 
by the (discretized) conformal time r, is 

N d (Ti)=V (n(Ti-i)-n(Ti)), (20) 

where n(r) = [^(t)]~ 3 is the number density of strings at 
time r. The energy- momentum tensor in Fourier space, 
0*, of these Ndin) segments is a sum 

%» = E . ( 21 ) 

771—1 

where 0™ is the Fourier transform of the energy- 
momentum of the m-th segment. If segments are un- 
correlated, then 

(e™0™> = ^777'(0™0™) (22) 

and 

(0™©™;'©"?") = W<W <ej£e£e#). (23) 

Here the angular brackets (...) denote the ensemble av- 
erage, which in our case means averaging over many re- 
alizations of the string network. If we are calculating 
power spectra, then the relevant quantities are the two- 
point functions of 0^„, namely 

N d (n) N d ( n ) 

KM P ) = ( E E e ^ ™/>- ( 24 ) 

m— 1 m' — l 
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Eq. 



allows us to write 



N d (Ti) 



(%u Q U) = E = Ndn){e%e%), (25) 



m— 1 



where 0* „ is of the energy-momentum of one of the seg- 
ments that decay by the time Tj. The last step in Eq. (Eq) 
is possible because the segments are statistically equiv- 
alent. Thus, if we only want to reproduce the correct 
power spectra in the limit of a large number of realiza- 
tions, we can replace the sum in Eq. ( pl| ) by 



(26) 



The total energy-momentum tensor of the network in 
Fourier space is a sum over the consolidated segments: 



K 



K 



(27) 



So, instead of summing over 2i=i-^d( r i) ~ l^ 12 seg- 
ments we now sum over only K segments, making K a 
parameter. 

For the three-point functions we extend the above pro- 
cedure. Instead of Eqs. (p3) and (p|) we now write 



N d (n) N d (n) N d (n) 

(%Mp e U) = ( E E E ©p©™'©;"/') 

m—l m' — l rn" — l 

N d (n) 

= E <<W?e#> = ^(n)(e« e^e; 1 ,) . (28) 

m— 1 

Therefore, for the purpose of calculation of three-point 
functions, the sum in Eq. @ should now be replaced by 

= [N d (n)} 1/3 %\ ■ (29) 

Both expressions in Eqs. ( p6| ) and (p9|), depend on 
the simulation volume, Vq, contained in the definition 
of Ndin) given in Eq. (20). This is to be expected and is 



consistent with our calculations, since this volume can- 
cels in expressions for observable quantities. 

Note also that the simulation model in its present form 
does not allow computation of CMB sky maps. This 
is because the method of finding the two- and three- 
point functions as we described involves "consolidated" 
quantities which do not correspond to the energy- 
momentum tensor of a real string network. These quan- 
tities are auxiliary and specially prepared to give the cor- 
rect two- or three-point functions after ensemble averag- 
ing. 



and a cosmological constant. Only the scalar contribu- 
tion to the anisotropy has been included. Vector and 
tensor contributions are known to be relatively insignif- 
icant for local cosmic strings and can safely be ignored 
in this model [ Q , M n. The plots are produced using a 
single realization of the string network by averaging over 

1 /3 

720 directions of k^. The comparison of G ; /j (or equiva- 

lently G^ 3 ) with its cosmic variance [cf. Eq. @] clearly 
shows that the bispectrum (as computed from our cosmic 
string model) lies hidden in the theoretical noise and is 
therefore undetectable for any given value of I. 

Let us note, however, that in its present stage our 
string code describes Brownian, wiggly long strings in 
spite of the fact that long strin gs a re very likely not 
Brownian on the smallest scales j22|. In addition, the 
presence of small string loops pqj and gravitational ra- 
diation into which they decay were not yet included in 
our model. These are important effects that could, in 
principle, change our predictions for the string-generated 
CMB bispectrum on very small angular scales. 
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FIG. 1. The CMB angular bispectrum in the 'diagonal' 
case (G^ 3 ) from wiggly cosmic strings in a spatially flat model 
with cosmological parameters ^CIDM; = ^baryon = 0-05, 
Qa = 0.65, and Hubble constant H = 0.65kms _1 Mpc _1 [up- 
per panel] . In the lower panel we show the ratio of the signal 
to theoretical noise \Ciu/(Tg l 1 ^ 3 for different multipole in- 
dices. Normalization follows from fitting the power spectrum 
to the BOOMERanG and MAXIMA data. 

The imprint of cosmic strings on the CMB is a com- 
bination of different effects. Prior to the time of recom- 
bination strings induce density and velocity fluctuations 
on the surrounding matter. During the period of last 
scattering these fluctuations are imprinted on the CMB 
through the Sachs- Wolfe effect: namely, temperature 



B. Results and discussion 



In Fig. we show the results for G^ 3 [cf. Eq. @]. It 
was calculated using the string model with 800 consoli- 
dated segments in a flat universe with cold dark matter 



"The contribution of vector and tensor modes is large in the 
case of global strings 
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fluctuations arise because relic photons encounter a grav- 
itational potential with spatially dependent depth. In ad- 
dition to the Sachs- Wolfe effect, moving long strings drag 
the surrounding plasma and produce velocity fields that 
cause temperature anisotropies due to Doppler shifts. 
While a string segment by itself is a highly non-Gaussian 
object, fluctuations induced by string segments before 
recombination are a superposition of effects of many ran- 
dom strings stirring the primordial plasma. These fluc- 
tuations are thus expected to be Gaussian as a result of 
the central limit theorem. 

As the universe becomes transparent, strings continue 
to leave their imprint on the CMB mainly due to the 
Kaiser-Stebbins effect . This effect results in line dis- 
continuities in the temperature field of photons passing 
on opposite sides of a moving long string.H However, this 
effect can result in non-Gaussian perturbations only on 
sufficiently small scales. This is because on scales larger 
than the characteristic inter-string separation at the time 
of the radiation-matter equality, the CMB temperature 
perturbations result from superposition of effects of many 
strings and are likely to be Gaussian. Avelino et al. H] 
applied several non-Gaussian tests to the perturbations 
seeded by cosmic strings. They found the density field 
distribution to be close to Gaussian on scales larger than 
1.5(Q.Mh 2 ) Mpc, where Qm is the fraction of cosmo- 
logical matter density in baryons and CDM combined. 
Scales this small correspond to the multipole index of or- 
der I ~ 10 4 . We have not attempted a calculation the 
CMB bispectrum on these scales because the linear ap- 
proximation is almost guaranteed to fail at such small 
scales, and because of increased computational cost for 
higher I multipoles. 

In summary, we have developed a numerical method 
to compute from first principles one of the clean- 
est non-Gaussian discriminators — the CMB angular 
bispectrum — in any active model of structure formation, 
such as cosmic defects, where the energy-momentum ten- 
sor is known or can be simulated. Our method does not 
use CMB sky maps and requires a moderate amount of 
computations. We applied this method to the compu- 
tation of some relevant components of the bispectrum 
produced from a model of cosmic strings and found that 
the non-zero non-Gaussian signal is unobservable even 
with forthcoming satellite-based CMB missions. Further 
computations and improvements using this method will 
be reported elsewhere. 



tin an extension of the Kaiser-Stebbins effect, Benabed and 
Bernardeau J5f| have recently considered the generation of a 
B-type polarization field out of E-type polarization, through 
gravitational lensing on a cosmic string. 
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APPENDIX A: PLOTTING THE BISPECTRUM 

We explain here our choice for the normalization of the 
angular bispectrum, as given in Eq. (|l^). The argument 
is an extension of the case of the angular power spectrum. 

We can express the two-point correlation function at 
zero lag in terms of the angular spectrum as follows 



while for the second sum we have 




(Al) 



In the small angular scale limit, the approximation that 
Ci is a smooth function of the multipole index I is well 
justified. We can then replace the sum by an integral and 
get 




l(2l + l)Ci. 



(A2) 



Now, dl/l — d(ln(l)) and therefore 1(21 + l)Ci is the con- 
tribution to the mean squared anisotropy of temperature 
fluctuations per unit logarithmic interval of I. In stan- 
dard practice, one usually plots 1(1 + l)Ci/2ir versus 
which for large I is proportional to 1(21 + \)Ci/Att. On 
small angular scales, then, this is oc l 2 C\. 

In the case of the three-point correlation function the 
situation is a bit more involved. Let us consider the skew- 
ness: 




Chhh 



(A3) 



We know that Ci^i^ is smooth in all three indices. We 
can split the skewness into three sums: the sum of terms 
where all li are different, the sum where only two of the 
three U are different, and the sum of terms where all U 
are equal. Omitting constant factors of 47r, we outline 
the same procedure as above for the two-point function. 
For the first sum we get 

) ' ( A4 ) 



ji/4i«li M (SSJ)(i 



1,1-2 



and for the third 
dl 



(A5) 



(A6) 



If one is interested in the diagonal terms Cm , then, fol- 
lowing the last equation, the relevant quantity to plot is 
given by 



1(21 + I) 3 / 2 [ I I I )(-',„. 





which is oc l 3 ^ 2 Cm at large I. 



(A7) 



7 



